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ABSTRACT 

From the rate of hydrogen ionization and the gamma ray flux, we derived the 
spectrum of relativistic and subrelativistic cosmic rays (CRs) nearby and inside 
the molecular cloud Sgr B2 near the Galactic Center (GC). We studied two cases 
of CR propagation in molecular clouds: free propagation and scattering of parti¬ 
cles by magnetic fluctuations excited by the neutral gas turbulence. We showed 
that in the latter case CR propagation inside the cloud can be described as diffu¬ 
sion with the coefficient ~ 3 x 10^^ cm^ s“^. For the case of hydrogen ionization 
by subrelativistic protons, we showed that their spectrum outside the cloud is 
quite hard with the spectral index <5 > — 1. The energy density of subrelativistic 
protons (> 50 eV cm“^) is one order of magnitude higher than that of relativis¬ 
tic CRs. These protons generate the 6.4 keV emission from Sgr B2, which was 
about 30% of the flux observed by Suzaku in 2013. Future observations for the 
period after 2013 may discover the background flux generated by subrelativistic 
CRs in Sgr B2. Alternatively hydrogen ionization of the molecular gas in Sgr 
B2 may be caused by high energy electrons. We showed that the spectrum of 
electron bremsstrahlung is harder than the observed continuum from Sgr B2, and 
in principle this X-ray component provided by electrons could be seen from the 
INTEGRAL data as a stationary high energy excess above the observed spectrum 


1. Introduction 


One of the important but still unresolved problems of astrophysics is spatial distribution 
of cosmic rays (CRs) in the Galaxy. We still do not have reliable understanding on how 
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CRs are distributed in different parts of the Galaxy and what is their spectrum there. In 
principle, this can be estimated from the characteristics of nonthermal emission generated 
by CRs in a range of wavelengths. The density of relativistic protons in the Galactic Disk 
(GD) can be derived from the observed intensity of gamma-rays, which are believed to be 
produced by proton-proton (pp) collisions while the density of subrelativistic CRs can be 


estimated from the emission of nuclear de-excitation lines and 
hydrogen (see, e.g.. 


Berezinskii et al. 


1990 


Ramatv et al. 


1979 


’ l-absorption lin es of ionized 


Oka et al. 


200511 . 


Several attempts were undertaken to derive the distribution of relativistic CRs in GD 


with energies abo ve 1 GeV from the gamma-ray data 




Tibaldo et al. 


'he distribution o 


CRs, derived 


Strong & Ma 

ttox 

1996 

and Fermi-LAT ( 

Abdo et al. 

2010; 

Ackermann et al. 

2011 


201311 measurements, showed that the density of CRs in the Galactic Center 


(GC) was higher than near Earth but their spatial distribution in the GD was much flatter 


Case & Bhattacharva 

199§; 

Green 

— S 

o 

CM 


20121 ) or pulsars (see, e.g.. 


Lorimer 


20M). 


A natural explanation could be an effective spatial mixing of CRs by particle scattering 
on magnetic fluctuations in the Galaxy if the GD is surrounded by a giant halo, in 


which CRs spenc 


Berezinskii et al. 


1988 


a SI 


gnihcant part of their lifetime before escaping from 


Bloemen et al. 


199011. H owever, numerical calculations showed (see, e.g.. 


he Galaxy (see 


Dogiel fc Urvson 


199311 . even in the most favorable case of a very extended halo, that 


this scattering (described as diffusion propagation) was unable to remove the signature of 
the source distribution. The CR distribution derived from the diffusion model was steeper 


than the distribution inferred from gamma-ray observations, alth ough it was f 


the source distribution. One of the explanations was suggested by 


Strong et al 


her than 


3 ) 


(120041 ) who 


showed that the problem can be solved by a variation in the Wco-^o-N{H 2 ) (metalicity) 
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scaling factor. Another explanation was snggested by iBreitschwerdt et al.l (120021) who 
assnmed that CRs left the disk faster from regions of a higher concentration of SNRs which 
will smoothen the CR distribntion in the GD in comparison with that of their sources. 

However, estimates of the CR density, derived from the gamma-ray data, are model 
dependent, because a simplihed model of the gas distribution in the GD is used in 
calculations. Besides, the adopted distribution of CR sources in the GD cannot be accepted 
as absolutely reliable because of the obscurity of dust and so on. 

Another possibility to dehne the density of CRs at different galactocentric radii (and 
even at different altitudes above the GD) is the analysis of gamma - ray emission from 


molecular c l ouds w ho se total mass i s know n, see 


Yang et al. 


(j2014aljb|); 


Tibaldo et al. 


Digel et al. 


(119961) : 


Hui et al. 


(120151) for the GeV energy range and 


( 120111 ): 


Aharonian et al. 


(120061) for the TeV energy range. A special case is the circumnuclear disk of molecular gas 


second and third Fermi LAT source catalog (see 

gamma-rav emission can be provided bv CRs g:€ 

Nolan et al. 2012: 

Abdo et al. 2015). This 

merated in the course of present or past 

activity of Sgr A* (see Chernvakova et al. 2011: 

Chernyshov et al. 

2014: Malvshev et al. 


20151) . However, analysis of the radiation from the source J1745.6-2858 is beyond the scope 


of the present work. We suppose to analyse this emission elsewhere. We note that 


van Eldik 


(120151) provided an excellent review on the origin of gamma ray emission from the GC. 


Recent analysis of gamma-ray emiss ion from lo c al mole cular cloud near Earth and 


the cloud Sgr B2 in the GC provided by 


Yang et al. 


(j2014aUbf) showed that the density of 


CRs was almost the same in these regions although the distance between them is about 8 
kpc. However, these estimations of CR density inside the clouds depend on how freely can 
energetic particles penetrate into molecular clouds. This question is one of the goals of this 
article. 
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The spectru m of CRs in th e energ y range below 1 GeV can also be estimated from the 
gamma-ray data. 


Dermer et al 


(120131) derived the spectrum of proton component of CRs 


in the local i nterstellar gas from the mid-latitude measurements of the Galactic gamma-ray 


emission see 


Abdo et ah 


20091) . These investigations showed that there was a flattening in 


the pr oton spectrum below several GeV. Similar conclusion was obtained by 


Neronov et al. 


(12012f) from investigation of gamma-ray emission generated in local molecular clouds. 


However, due to the energy threshold of pp reaction, the gamma-ray data that is useful for 
deriving the spectrum of CR protons are those with energies above hundred MeV. 


Bouchet et al.l (120111) derived the spectrum of Galactic electrons in the range between 
1 and 5 GeV from the intensity of diffuse hard X-ray emission measured by INTEGRAL. 
They showed that a signihcant fraction of the Ridge emission in the range above 50 keV is 
produced by the inverse Compton scattering of these electrons on background photons. 


Almost forty years ago iRamatv et al.l (119791) suggested to measure a flux of nuclear 
de-excitation lines with ~ 0.1 — 10 MeV generated by subrelativistic CRs. In 
order to estimate their density in the range between several MeV and hundred MeV. 


However, the sensitivity of gamma-ray telescopes was not (and sti 
to detect these lines in the diffu se gamma-ray spectrum (see e.g. 


1 is not) hig 


Benhabiles-Mezhond et al. 


20131) 


Dogiel et al. 


1 enou gh 


2009b 


Another very interestin, 
GC region was obtained by 


; results about the par ameters of low energy cosmic rays in the 


Nobnkawa et al. 


(I 2 OI 5 I ). They found an excess of the 6.4 keV 


line emission in the near east of the GC. They concluded that this excess is due to iron 
atom bombardment by protons with energies 0.1 - 1000 MeV. The estimated energy density 
of these protons is about 80 eV cm“^ which is almost two orders of magnitude higher than 
the energy density of relativistic CRs near Earth. 


In principle, the spectrum of low energy CRs in the range below 100 MeV can also 
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be derived from the ionization rate of hydrogen in the interstellar medium. The point 


is that subrelativistic electrons and protons ionize hydrogen molecule s effective^ 
ionization rate can be d erived from the IRabsorp 


]ion lines (see 


Indriolo fc McCall 


2012 


see also the review of 


Dalgarno 


Oka et al 


2001 


re 


(120061) 1 that allows one to estimate 


the dens ity of subrelativistic CRs insid e molecular clouds and in the intercloud medium 


(see e.g. 


Dogiel et ah 


2002 


2013 


20i4 . 


The average ionization rate in the 100 pc central region of the Galaxy, measured by 


Oka et ah 


(120051) . is about (3 — 5) x 10 s while this value for the cloud Sgr B2 is 
depleted and equals 4 x 10“^® s“^ (iVan der Tak et al.ll2006l) . This means that penetration 


of subrelativistic CRs into Sgr B2 is not complete, and the average density of these CRs in 
the intracloud region is higher than inside Sgr B2. Similar effect is observed for the clouds 


-16 ^-1 


s ^ in the diffuse c 


outside the GC where the ionization rate is about 3.5 x 1C 
while this value is 3 x 10“^^ s“^ in dense molecular clouds (|Indriolo fc McCall 


2012 


ouds, 

)• 


The spectrum of subrelativistic CRs inside molecular clouds depends on processes 
of particle propagation there. Several models are used in order to estimate CR density 
in the clouds. One of them is extrapolation of the model of CR propagation in the 


Dogiel et al. 

20 Q2; 

Cabici et al. 

o 

o 

CM 


Protheroe et al. 


2008). Alternatively, 


Padovani et al. 


(l2009h derived the spectrum of subrelativistic CRs in the local interstellar medium from 
the rate of hydrogen ionization in the local molecular clouds assuming that they propagate 


freely (without scattering) inside 
may be essential in this case (see 


he clouds althoug. 


Padovani fc Galli 


1 the effect of large scale magnetic held 


20111) . We discussed processes of CR 


penetration into the clouds and show that this process may differ strongly than supposed 
in previous publications. We try to estimate also the spectrum of CRs in the GC from 
the observed gamma-ray hux and the ionization rate of hydrogen in Sgr B2. Propagation 
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of CRs in partially ionized medium is a rather complicated process, which includes CR 
ionization of the medium, hydromagnetic waves generation and damping, density and 
magnetic structures of the medium. We address the density structure and ionization state 
of Sgr B2 in Section 2. Section 3 discusses excitation of hydromagnetic waves by CRs, wave 
damping by ion-neutral collisions and diffusion of CRs in the cloud, and Section 4 considers 
the magnetic structure of the cloud. Base on these, we calculate the CR spectrum for two 
models of the cloud in Section 5. Section 6 gives some consequences of the model. Section 
7 discusses the possibility of ionization by electrons. Finally, in Section 8 we provide a 
summary. 


2. Hydrogen parameters of Sgr B2 and energy losses of CRs 


Parameters of CR propagation depends on the medium parameters. The strength 


of magnetic held 

0.5 mC (see, e.e.. 

n Sgr B2 was estimated from the Zeeman splitting; and was about 

Cm 

:cher et al. 

1996 

). The cloud Sgr B2 is a very massive complex. 

Aharonian et al. 

( 

2006 

) estimated its total mass as (6 — 15) x 10® Mq for a 0.5° x 0.5° 


(75 X 75 pc^) region surrounding Sgr B2. The gas distribution in the cloud Sgr B2 is highly 
uncertain. 


Lis fc GoldsmithI ( 1989 


1990, 


1991 1 and 


Goldsmith et al. 


(119901) derived a two- 


component density distribution: an envelope with a constant H 2 density about 1800-3500 
cm“^ extending to the outer radius about 22.5 pc and a central core with the central density 
4 — 9 X lO'^ cm“^. According to their analysis the density distribution can be approximated 
as 

, , , ni+712, if r < To, 

nnAr) = 


ni{^y + n2, if R>r>ro, 


( 1 ) 

























with ro = 1.25 pc, R = 22.5 pc, ni = 5.5 x 10^ cm“^, 77-2 = 2.2 x 10^ cm“^, and a = 2. The 
H 2 column density at the center is, N{H 2 ) = 2.6 x cm“^, and the total mass of Sgr B2 
they took as 6.3 x 10® Mq. 


Protheroe et al.l (120081) expressed the gas distribution in Sgr B2 in the form 


M, 


nH2 — 


H 2 


2mH (\/^i 


TTCT 


■exp(-r^/2(T^), 


( 2 ) 


where Mh 2 is the total mass of Sgr B2 and a = 2.75 pc. The density is about 10® cm“^ at 
the center of the Sgr B2 complex, and decreases to 10 cm“® at a radius of ~ 12 pc, which 
they considered to be its outer radius. The central column density is about 2.5 x 10 ^ 
and the estimated mass of Sgr B2 is about Mh^ ~ 2 x 10® Mq (see 


is about 2.5 x 1( 

l^^cm' 

Protheroe et ah 

2008) 


Below we analyze processes of hydrogen ionization and gamma-ray emission for these 
two extreme distributions of hydrogen in Sgr B2. 

CRs lose effectively their energy inside the dense molecular clouds. In the non¬ 


relativist ic ener gy range the rate of losses is determined by ionization (see 


Ginzburg 


19891 ) 


Havakawa 


19641 : 


(dE\ 


2'Ke^nH , f mlc^W„ 

In 


(3) 


meCj3{E) \ A7re'^h?n 

where rin is the density of background gas, me is electron mass, ITmax is the highest energy 
transmitted to an ambient electron, and P{E) = u/c, where /5 is the particle velocity u in 
the unit of light speed c. The characteristic time of the losses is 

f dE 


n = 


{dE/dt)i 


(4) 


In the relativistic energy range protons lose their energy by collisions with the ambient 
gas {pp collisions). The characteristic time of the process is 


'^pp — ^'^H^pp 1 


( 5 ) 
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where the cross-section of pp-reaction can be found in e.g. 


Kamae et ah 


fl2nnfih 


One of the key point for the problem of CR penetration into dense molecular clouds is 
how they propagate through the gas. This problem is discussed in the next section. 


CR diffusion inside molecular clouds 


In the standard m odel of CR propagation through the interstellar medium (see, e.g., 


Berezinskiietah 


19901 ) the spatial diffusion coefficient can be presented in the form 

2 


D 


u 


v(E)' 


where ^{E) is the frequency of particle scattering by resonant MHD-waves: 

6H{kf 


V ~ UJii- 


-^0 


( 6 ) 


(7) 


Here 6H{k) is the strength of a magnetic fluctuation with the wave number /c, Hq is the 
large scale magnetic held, ojh = eHo/mc'j, 7 and m are the gamma-factor and the mass of 
a scattered particle, respectively. The wave-particle interaction is resonant. A particle with 
the energy E is scattered by waves whose wavelength is about the particle Larmor radius, 
A = 27r//c ~ tl^E). The amplitude of magnetic huctuations, hif, is supposed to be much 
smaller than the strength of the large scale magnetic held EIq in the interstellar medium. 


5H{k) < Ro • 


( 8 ) 


For the wave spectrum of huctuations W{k) 0 : k where 5H{kY = kWjk), the spatial 


dihusion coefficient D of relativistic CRs is energy dependent as (see 


Berezinskii et al 


19901) 


D{E) oc E 


2—k 


(9) 


This model of CR dihusion in the interstellar medium is often extrapolated to the case 
of CR propagation inside molecular clouds. However, the process of particle propagation 
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nearby and inside almost neutral and dense gas of molecular clouds is quite different in 
comparison with the standard process of CR diffusion in the interstellar medium. There are 
several specihc mechanisms which determine this process: 


1. Dense molecular clouds absorb CRs with energies whose lifetime Tit{E) (determined 
by Equation (j4)) or ([5])) is smaller than the time of particle propagation through 
the cloud, Tpr- For ballistic propa gation in the cloud Tpr = R/n and for diffusion 


propagation Tpr = /Dc (see e.g. 


Morhll 


1982 


Dogiel fc Sharov 


19851 ). Here R is a 


cloud radius, u is the particle velocity, and Df. is the coefficient of spatial diffusion in 
the cloud. 


For energies determined by the condition, Tit{E) < Tpr{E), a flux of CRs to the 
cloud arises be cause of CR absorption. This CR flux generates MHD-waves with the 


increment (see 


Kulsrud fc Pearce 


196911 


7(k.r) = ^ Jdpud - , (10) 


P 

Here H is the magnetic held strength near the cloud, va = H/yjAnpi is the Alfven 
velocity, pi is the density of ionized gas, /(p, r,/i) is the CR particle distribution 
function, p is the particle momentum, p is the particle pitch-angle, k is the wave 
number of a magnetic huctuation excited by CRs, and = H • k/|H|. 


2. On the other hand, ion-neutral friction in the dense low ionised molecular clouds 
damps MHD-waves with frequencies oj = fcVh < Pin, where Va = HU/AnPi is the 


Alfve nic velocity, pi is the density of ionised component (see e.g. 


Kulsrud fc Pearce 


19691 ). and the decrement of wave absorption is 

rrin 




{rrii + rrin) 


nnio'inVH) ■ 


( 11 ) 


Here ain is the cross-section of ion-neutral collisions, is the mass of neutral 
particles, rrii is the mass of ionised particles, and vh is the thermal velocity of 
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hydro gen. For the condition of molecular clouds {av) = 2-10 ®cm^s ^ fsee lPinto et ah 


200811 ■ 


Dogiel fc SharovI (1198511 provided detailed analysis of CR propagation nearby 


molecular clouds and showed that MHD waves are excited far away from the clouds, 
but they are completely damped inside the clouds. This means that CRs should freely 
propagate inside the clouds without scattering (but with energy loss). We notice also 


that similar analysis of CR propagation nearby the c 


there was provided latter by 


ouds and MHD-wave excitation 


Everett fc Zweibell ( 2011 1 and iMorlino fc Gabicil ( 2015 1. 


3. However, as obser vations showed, the neutral ga s in the clouds is highly turbulent 


(see the review of 


Hennebelle fc Falgarone 


2 OI 2 I) . This turbulence of neutral gas 


excites forcec 


(IDogiel et al 


magn etic fluctuations through interaction with the ionized component 


198711 . These fluctuations prevent free particle propagation in the 


clouds. This process is discussed in the next section. 


4. Magnetic field structure inside molecular clouds 


by 


Detailed analysis of ma gnetic held structure in turbulent molecular clouds was provided 


Istomin fc Kiselev! (120131) . Below we present a brief review of their analysis and apply its 


results to the case of Sgr B2. For an homogeneous and isotropic medium they considered 
the correlation function of the velocity held (nj(x, f)nj(x + r, f)) as delta-correlated in time, 

(nj(x,t)nj(x-F r,f)) = Vij{r)Tc6{t - f), (12) 

where Tc is the characteristic time of hydrodynamic turbulence of the neutral gas. 

From the conditions of the tensor symmetry and for the incompressible liquid (V-v = 0) 
the correlation tensor Vij{r) can be presented in the form 


-'o 


(r) = 2V{r)5ij + r 


dV{r) 

dr 




TiTj 




), 


(13) 
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where V (r) is an unknown function, which will be derived below from the observed spectrum 
of turbulence in molecular clouds. 


In this medium the correlation tensor of magnetic fluctuations bij similar to the velocity 
correlation tensor, but the average is taken at the same time moments, 

bij{r, t) = {bi{x, t)bj{x + r, t)) = 2Q6ij + (^5ij - ^ j , (14) 

where the function Q is a function of r and t in the general case. The amplitude of magnetic 
fluctuations bo is dehned as 6 q = ( 6 j(x) 6 j(x)). 


If there is a large scale magnetic held Ho in the medium, then t he equation for the 


corre lation function Q can be derived from the MHD equations (see 


Istomin fc Kiselev 


201311 


1 dQ{r) 

2Tr dt 


1/(0) - V{r) + 


dVdQ 1 fdV 

-H-4-^ r - 

dr dr r \ dr dr'^ 


^ I 
0 


d^Q 4 dQ 
dr‘i r dr 


(15) 


where pi is the density of ionised component of the gas. 


The spectrum of turbulent velocities in molecular clouds was derived from measurements 
of the CO and NH 3 line Doppler broadening. The turbulence h as a power-law Kolmogorov- 


l ike spectrum in a very broad range of scal es from supersonic flLarsqn 


(IMvers 


19831) regions (see also the review of iHennebelle fc Falgaronell2012h : 


198 ill to subsonic 


v{L) = 1.1 T"(pc) km s ^ where a ~ 0.3 — 0.5 , (16) 


where v is the velocity of turbulent motions and L is its scale (0.01 < L < 300 pc). Then 
from Equation flT^ the correlation function of velocity V can be derived. With the known 
V we can dehne the correlation function of magnetic huctuations Q from Eq. flT^ . Below 
we obtain this function for the stationary case when dQ/dt = 0. 
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Here we use dimensionless units for the magnetic field strength; 

H= ^ H. 

Here Uq is the the velocity of turbulent motions in the scale L < R. 


(17) 


If dimensionless large scale magnetic held Hq « 1, the energy of magnetic held 
huctuations is concentrated near the correlation length, Lcorri which is much smaller than 
the cloud size R. Then the a mplitude of magnetic huctuations, 5H, is much larger than Hq 


see 


Istomin &: Kiselev 


20131), 


5H^ Ho, 


( 18 ) 


(cf. Equation ([8])). 


By dehnition the correlation length is 


Lcorr - j (^i(0)&i(r))dr , 
° 0 


(19) 


and can be estimated if the function Q is known from Equation (|T5l) . 
Then from numerical calculations of Equation flT^ we obtain that 


r/R = 0.9-ffo ) 


SH = 


( 20 ) 

( 21 ) 


We apply this theory to the molecular cloud Sgr B2. For estimates we use the following 
parameters of the cloud and the intercloud medium. The gas density in the cloud uh — 10 ^ 


cm_f, the ionization degree in the dense molecular clouds njj+/nH 2 ~ 10 ® (s ee e.g. 


2nn6h, the radius of the Sgr B2 region emitting gamma-rays is i? = 7 pc 


19961: 


2014b), the total magnetic held stren . 


Crutcher 


1999 


Crutcher et al. 


GC is about TTo ~ 50 — 100 pG (jCrocker et ah 


;th in Sgr B2, Hq -|- 6H ~ 550 pG f Crutcher et al. 


see 


Oka 


Yang et al 


2010h. the avera ge large scale magnetic held in the 


2010h . For the core Sgr B2, where the 
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gas density is uh ~ 10^ cm and SH/Hq ~ 5.5, Equation fl20|) gives Lcorr — 0.4 pc, and 
Equation fl^T]) gi ves Vn = 7. 3 km s ~^. This estimates of the turbulent velocity is close to 


Crutcher 


fll999l) for Sgr B2. In the outer envelope with ~ 10^ — 10^ 


that presented in 
cm“^, we expect that hif Hq- 

In a magnetic held H = 0.55 mG the Larmor radius, of particles with energies 
E < 10^ GeV is smaller than Lcorr- In this medium propagation of magnetized particles 
along t angled magnetic h eld lines, can be described as dihusion with the coefficient (see for 


details 


Do del et ah 


1987h 


Dc{E) 


cP{E)Lc 


cLr 


{E/mpC^Y + 2{E/mpC^) 


2 ^ \l {E/mpC^f + 2{E/mpC^) + 1 ’ 

where /3{E) = u/c and u is the particle velocity . It follows from Equation fl22|) that 


( 22 ) 


Dc{E) oc 


\/E, 


if u < c, 


constant, if u c, 


(23) 


(cf. Equation ([9])). 


5. Spectrum of CRs inside and outside Sgr B2 



ionization rate of hydrogen, is 


c= aHuN{E)dE, 


(24) 


where an is the ionization cross-s ection of the molecular hydrogen by proton impact (see 


Rudd et ah 


1983 : 


Tatischeh 


20031) . and u is the velocity of CR particles and N{E) is their 


spectrum. 
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As we mentioned in Section [T] the ionization rate in the diffuse molecular gas, Co is 
about 3 X 10~^^ s~^, while inside t he co re of Sgr B2 the rate, C is o ne order of magnitude 


smaller. 


Dogiel et al 


(12011 


2013 


20141) and 


Tatischeff et al 


(120121) presented arguments 


in favour of hydrogen ionization by subrelativistic elec trons, while the al 


Yusef-Zadeh et al 


ernati ve process 


( 2013 ). Below we 


of ionization by high energy electrons was discussed in 
discuss these models. 

For calculations we approximate the gas distribution in Sgr B2 according to 
Equation ([I]) or (E]): 


Mode l I; Sgr B2 consists of the two components as it was derived by Ibis fc Goldsmith 
fll99lf) . a dense core with density ~ 10^ cm“^ and radius tq = 5 pc surrounded by 
a diffuse component with density Ud ~ 10^ — 10^ cm“^ and radius i? = 20 pc; 


Model II: The mass o 


Protheroe et al. 


the cloud is concentrated in a dense core (^as it follows from 


(120081) ) with density Uc ~ 10^ cm ^ and radius tq = 5 pc without an 


outer envelope. 


For both distributions the total mass of the cloud is between (2 — 6) x 10® Mq and the 
gas column density in the direction of the cloud center is ~ 10^“^ cm“^ as follows from 
observations. 

From Section 0] it is clear that magnetic fluctuations with 6H > Hq are excited in the 
dense core where the kinetic energy of turbulent motions is high enough. Only there we 
expect diffusion propagation of CRs with the coefficient (123|) . In the envelope, particles 
propagate without signihcant scattering because the fluctuations are damped by ion-neutral 
friction. 

We assume that the spectrum of subrelativistic protons, N{E), in the intercloud 
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medium can be presented as power-law, 

- E), (25) 

where E is the kinetic energy of subrelativistic protons, -Emax is the maximum energy of 


N{E) = K e{E„ 



The magnetic helds are nonuniform inside molecular clouds and their strength 
correlates with the gas density as 


H (X , 


(26) 


see 


Crutcher et ah 


(l2010l) . We cannot say whether it is due to the frozen effect of the 


magnetic helds into the gas only or this may be due to more effective generation of magnetic 
huctuations in dense cores by the neutral gas turbulence as we assumed above. 


Padovani fc Gallil (120111) assumed a convergent structure of magnetic held line nearby 


molecular clouds and analysed the ehect of CR mirroring. They showed that this ehect 
reduced the CR density in the clouds and the ionization rate there by a factor of ~ 4. 

In the following we neglect variations of the large scale magnetic held in the outer 
envelope. We assume that inside the outer envelope subrelativistic CRs propagate through 
the medium without signihcant scattering, i.e., along held lines of the large scale magnetic 
held Ho with their own velocity. This is a big simplihcation of the process, but it provides 
an upper limit for the density of CRs penetrating into the envelope. In this work we 
consider stationary distribution of CRs in molecular clouds. Unlike the case of the source 
J1745.6—2858 near Sgr A* (right at the GC), we do not have evidences indicating a time 


varying CR source at or near t 
for the time-dependent case in 


le location of Sgr B2. N evertheless, one can hnd a solution 


Chernyshov et ah 


(120141 ). 
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We describe CR propagation inside the onter envelope as 

= 0 , 


d \fdE\ 

ae [{lu 


(27) 


with the spectrnm of CRs at the outer boundary of the envelope (x = 0) as presented by 
fl2^ . The goal is to calculate the spectrum of subrelativistic CRs that reach the surface of 
the core at ro = 5 pc. 


The solution of Eq. fE7|) is 


N(x,E)^K{ — 


1/2 


(- 

\E„ 


ax\/m 

-^0 


^+0.5 

2 


e{E, 


2 

max 


E'‘ - . 


(28) 


Here we presented the term of energy losses in Eq.(l3]) as dE/dt = a/\/E, and the mean free 
path of a subrelativistic particle with mass m and energy E is x E / [dE / dt)i^j2E/m = 
E'^/a^/m. 


The characteristic lifetime of a subrelativistic proton with energy E is 


dE 


m, 


,^ 3/2 


r = 


{dE/dt)i ^ 


(29) 


The energy loss gas column density, which a subrelativitic proton with energy E passes 
through during its lifetime, is 


Lh — - 


1 


3 m vre'^A 


(30) 


In Figure [T] we show the energy loss gas colu mn density 


derived from equations o 


Mannheim & Schlickeisei 

(1 

99i: 

1 ). F 

taken from 

Lis & Goldsmit’ 

3 C 



Hayakawal fll96 


or electrons anc 


3; 


protons 


GinzbnrgI fll989l ) and 


1990, 


gas column density of the outer envelope 


19911) . which is presented by Equation ([T]), we 


can estimate the energy of protons which can reach the surface of the central core. For 
Lh — 10^^ cm“^ this value is about E > 20 MeV . 
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Fig. 1.— Energy loss gas column density, L^, passed by a proton (solid line) and an electron 
(dashed line) with energy E. 


As an example, in Figure [2] we show the expected spectrum of protons at the boundary 
of the central core Nc{E) for Lh = 10^^ cm“^ if the spectrum of protons in the intercloud 
medium is N{E) = K{E/Eq)~^-^. 


According to Section HI inside th e core CR propagation is described as diffusion with 


the coefficient (see 


Dogiel et ah 


1987h 


Dr = 


uLr 


(31) 


where L^orr = 0.4 pc. 


The distribution function of protons, Nc{E), is described by the equation 


r 2 y dr ) dE [ \ dt / ^ 


0 , 


(32) 


where {dE/dt)i is the rate of ionization losses in the dense core, calculated for the density 
uh ~ 10^ cm“^. 















- 19 - 



Fig. 2.— Expected spectrum of proton at the boundary of the central core Nc{E) for Lh = 
10^^ cm“^ and the spectrum of protons in the intercloud medium N{E) = K{E/Eq)~^'^. 


The boundary condition on the core surface, r = tq 

_ j+0.5 

E \ 1/2 " 


^c{E) = 


E\ a{R — rQ)^/m 

~Eo) ^ ^ 


(33) 


where R and ro are the radius of the outer envelope and the central core (see Model I). We 
can ignore the unknown constant K from the ratio of ionization rates outside and inside 
the core. 


The ionization rate averaged over the scale Lh of the outer envelope is (^ = 3 x 10 


-15 


s \ 

J dr J (jHuNc{E,r)dE . 

R 

The ionization rate inside the central part of core is Cc = 3 x 10“^® s“^. Then 

R Emax. 

Jdx J dEN{x,E)ai{E)u{E) 

C '1^0 ro Eniin 

_ _ 

Cn (R — rn) ^0 Emax 

Jdr J dENc{r,E)ai{E)u{E) 

0 ^min 


(34) 


( 36 ) 
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From Equation fl35l) we can derive the spectral index 5 as a function of E mav for a given 
value of Lh in the outer envelope, where = uh^R — no), i? = 20 pc and tq = 5 pc. The 
constant K is calculated then from the ionization rates inside or outside the cloud. 

In Figure |3] we show the spectral index 5 for different values of nn- The case uh = 0 
corresponds to the case where the whole mass of Sgr B2 is concentrated in the core (i.e. 
Model II). One can see from the hgure that for very hard spectra with 5 = 0.2 almost all 



Fig. 3.— The required maximum energy of protons Emax as a function of the spectral index 
5 and the gas density in the onter envelope, when nu 7 ^ 0 in the outer envelope (Model I), 
and for = 0 in the outer envelope (Model II). 

particles reach the core surface. The required maximum energy E m ax of protons is in this 
case about 0.2 GeV independent of the gas density. 

The derived spectrum of subrelativistic protons is shown in Figure |4] for different values 
of Uh in the envelope and i? — tq = 15 pc. The derived values of E m ax , the density of 
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hydrogen uh, and the energy density of snbrelativistic protons needed for the ionization 
rate in Sgr B2 are also shown in the hgure. 


The relativistic component of CR spectrnm inside Sg r B2 can be derived from the 


observed flux of gamma-rays there (see 


Yang et ah 


2014bl) . Gamma-ray flux generated by 


relativistic protons in molecular clouds can be calculated by 

R 

da 


F,{E,)=A7rc I nnirydr / N{E,r)^{E, E,)dE , 


where nnij) is the hydrogen distribution, da/dE^( E, E^) is the c 


(36) 


ifferential cross-section for 


20061) . and N{E,r) is the 


gamma-ray production in proton-proton collisions (IKamae et ah 
distribution of relativistic protons in Sgr B2. 

If the mass of Sgr B2 is concentrated in the envelope (i.e. Model I), then according to 
our analysis relativistic protons hll uniformly this region and their density is the same as 
in the intercloud medium. It is clear from Figure [2] that even protons with E ~ 0.2 GeV 
can penetrate through the outer envelope. Protons with energies E >1 GeV pass through 
the envelope witho ut absorption {ji^Tpp > Tpr). The spectrum derived for this case by 


Yang et al. 


(j2014alJbl ) is shown as shaded areas in Figure 01 Variations of the spectrum are 
governed by the uncertainties of the mass of Sgr B2, which is in the range of 2 — 10 x 10® 


Mq. The energy density of relativistic protons is then within 1 to 4 eV cm 

It is interesting to notice that in this case the required energy density of subrelativistic 
protons is more than one order of magnitude higher than that of relativistic protons. 

If the mass of Sgr B2 is mainly concentrated in the core (i.e. Model II), the distribution 
of relativistic protons in the cloud is calculated similarly to Equation 


Tppiy 


= 0 , 

or \ or 


(37) 


where ~ cLcarrl‘^ =constant and = {nucapp^E)) ^ is the characteristic time of pp 


















collisions. The boundary condition is 


Nc{E) = N{E), 


(38) 


where N{E) is the spectrum of relativistic protons in the intercloud medium, which is 
supposed to be a power-law. 

In this case protons are distributed in the core nonuniformly (rpp < Tpr) and their 
spectrum can be derived from Equations (15B]) & fl57|) . Results of the calculations are shown 
as the solid line in Figure 01 One can see that the energy density of relativistic protons is 
about one order of magnitude higher than that near Earth, but it is still smaller than that 
of subrelativistic CRs. 

From Figure m it is reasonable to conclude that relativistic and subrelativistic protons 
have different sources in the GC because their spectra have a gap in intensities in the 
range about several hundred MeV and they do not match smoothly with each other at 
these energies. An exception is the case where the whole mass of Sgr B2 is concentrated 
in the core (the solid line in the subrelativistic energy range and the straight line in the 
relativistic energy range). It seems from the hgure that in this case the distribution of 
subrelativistic and relativistic protons can be described by a unihed spectrum with a break 
of the spectral index from 5 ~ —3 in the relativistic energy to 5 ~ 0.2 in the subrelativistic 
energies. Similar spectrum can be formed in the intercloud space, e.g., by ionization 
losses (see Equation ([3])) if the spectrum of protons injected by sources has a cut-off, 
Q{E) = KE^6{E — Emin). For the equation 



(39) 


the solution for E < E^^ is 


N{E) oc . 


( 40 ) 


For more details of this solution see iDogiel et al.l fl2009bl) . 
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6. By-products of the model of hydrogen ionization by subrelativistic protons 

With the spectrum of protons derived in the previous section we can analyse whether 
this model can be proved or disproved from observations. 


6.1. Emission of the 6.4 keV iron line 


First we would like to mention that subrelativistic protons not only ionize hydrogen 
molecules in the clouds but also atoms of heavier elements, e.g., iron, whose ionization 
of K-a shell is observed as the 6.4 keV line emission. It was proved that thi s emission is 
gener ated by a front of hard X-ray photons emitted in the past by Sgr A* (see 


Kovama et ah 


1996 ) whose luminosity was much higher in the past t 
20131 ). The emission shows time variability (see, e.g.. 


ran at present 


Kovama et al. 


(see e.g. 


Rvu et al. 


20081 ), and it had 


decreased by one orde r of m agnitude from the peak period in 2000 to 2013 fllnui et al. 


Nobukawa et al. 


2011 


2009 


2014J ). It is expected that for about twenty years of observations of 


Sgr B2, the front of X-rays have to leave (almost) the cloud. Then we may expect that a 
background flux of the li ne emission from S gr B2 generated by CRs can be seen at present 


or in the near future (see 


Dogiel et al. 


20111 ) . 


Temporal variations of the 6.4 keV Sgr B2 flux in the period 2005-2013 are presented 
in Table [U 

Below we compare these Suzaku data with model estimations. From the spectrum 
shown in Figure 0] we can estimate the expected 6.4 keV line emission generated by the 
protons: 

1 


R 


h.i — 


d2 

^GC 


nH{j')r^dr / N{E,r)vaQAdE, 


( 41 ) 


0 E 

where the cross-section of 6.4 keV photon production by proton impact is from 
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Table 1: The 6.4 keV flux from Sgr B2 region of 2.0' radius. 


Year 

Fel-KcK®'’^ 

2005 

(18.0 ±1.3) X lO-'^ 

2009 

(7.4 ±0.8) X 10-5 

2013 

(3.0 ±0.9) X 10-5 


“Errors at 90% confidence levels. Background is taken into account. 
''Observed flux in the unit of photons cni“^ s"'^. 


Note. — The data were derived from 


Nobukawa et al. 


( 2011 


201J). 


Tatischefll ( 20031) . 


The expected intensity depends weakly on model parameters and is about 


Iqa ~ (3 — 5) X 10 ph cm ^ s ^, 


(42) 


where rj is the abundance of iron atoms with respec 
?7 in the Sgr B2 is in the range 1.3 — 1.9 (see, e.g., 


) to solar. The es 

imated iron abundance 

Revnivtsev et al. 

2004; 

Nobukawa et al. 


20111) . The estimated flux of the 6.4 keV line in Eq. 


is three times smaller than it was 


in 2013. Hence in the framework of the model, subrelativistic protons might produce in 
2013 about 30% of the total flux from Sgr B2. We cannot exclude that at present, the front 
of hard X-rays has left Sgr B2 and this emission is generated by subrelativistic CRs. Then 
it should drop down to the level predicted by Eq. fl42|) . Future analysis of the X-ray data 
for the period after 2013 may confirm our conclusion. Restrictions on processe s of 6.4 keV 


line p roduction can be obtained from measurements of the line broadening (see 


Do del et al. 


19981) . We hope that this crucial information to constrain the origin of the 6.4 keV emission 


from Sgr B2 will be obtained by ASTRO-H (see 


Koyama et al. 


20M). 


























25 


6.2. Hard X-ray continuum 

Hard X-ray photons and subrelativistic CRs also produce a flux of continuum emission 
by Thomson scattering of photons ( Compton echo) or by bremsstrahlung losses of charged 


particles. Observations presented by 


Terrier et ah 


fl2010l) showed time variability of the 


X-ray flux from Sgr B2 in the range from 20 to 100 keV which correlated nicely with the 
variability of the 6.4 keV line emission from Sgr B2. This confirmed that their common 
origin is due to X-ray irradiation from Sgr A*. 


The X-ray emission produced by inverse bremsstrahlung of protons can be calculated 


from 

R 

h{.E^) = —^ f nH{r)r^dr j N{E,r)^^{E, E^)dE , 

^GC J J ^ 

0 E 

where a^/dE^ is the cross-section of inverse bremsstrahlung 


(43) 


da. 


dE. 




3 he KrUeC^ J EE. 


— \/ E — E^ 


E^ 


(44) 


where E = 


from ([Hang 


m^Jmp)E for protons and E = E for electrons. Below we adopt the cross-section 


19971 ). Subrelativistic p rotons wit 


with energies E^ < {me/mp)E (see 


1 the energy E emit bremsstrahlung photons 


Havakawa 


19641) . Therefore for the cutoff energy of 


subrelativistic protons about 0.2 — 0.5 GeV we expect that they generate hard X-ray 
emission in the range T'a; 10 — 30 keV. 

The calculated X-ray fluxes depending of the Sgr B2 mass are about 


= (4 — 6) X 10 erg s ^ cm ^ for 2 keV < E^ keV, 


= 7.5 X 10 erg s ^ 


cm 


-2 


for 20 keV < E^ < QQ keV. 


(45) 

(46) 


One can see that the model e stima tions are significantly below the experimental data 
presented in 


Nobukawa et ah 

(2011 

) for the range 2 — 10 keV and in 

Terrier et ah 

(2010) 
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for the range 20 — 60 keV. However, the continuum X-ray flux shows strong decrease in 
the Suzaku (2-10 keV) and the INTEGRAL (20-100 keV) ranges. For the period from 


2003 to 2009 the flux dropped down by a factor 0.4 flNobukawa et al 


2011 


Terrier et al 


20101) ■ Thus, one can expect that the continuum emission produced by subrelativistic CRs 
(Equations fj45l) & fj46l) ) may be seen in near future. 


6.3. De-excitation nuclear lines 


Subreletivistic protons and nuclei with energies between a few MeV and several 
hundred MeV generates in the interstellar medium e mission of nucle ar gamma-ray lines 


in the range 

Tom 0.1 to 

10 MeV fsee the review of 

Tatischefl 

2003 

). Estimations of 

Dogiel et al. 

I2nn9b 

) and 

Benhabiles-Mezhoud et al. 

(2013 

) showed, however, that this line 


emission cannot be detected even in the GC because the sensitivity of existing gamma-ray 
telescopes is not high enough. It seems nevertheless that the discovery of the gamma-ray 
lines will be possible with new generation telescopes. Bearing this in mind, we estimated 
the flux of G and O lines from Sgr B2, expected for the spectrum of subrelativistic protons, 
derived in Section |5l For the double solar abundance in the GG the flux of the O line from 
there is about 1.38 x 10“^ ph cm“^ s“^ and the G line 8.34 x 10“® ph cm“^s“^. 


6.4. Secondary electrons and radio emission from Sgr B2 


Relativistic protons generate from pp collision not only gamma-rays, as discussed in 
Section 0 but also secondary electrons, which should produce a synchrotron radio emission 
in the relatively strong magnetic helds of Sgr B2. From the observed gamma-ray flux we 
can estimate the density of secondary electrons and the flux of radio e mission using the 
equations for these processes presented, e.g., in 


Berezinskii et al. 


(Il99nh . The flux of radio 
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emissi on from Sgr B2 in the range 330 MHz to 2368 MHz was presented in Hones et al. 


( 120111 ) who concluded that they did not hnd any evidence for non-thermal emission from 


Sgr B2 and the observed radio flux from there is of the thermal origin. 

We estimated the synchrotron component generated by secondary electrons. As an 
upper limit for the density of secondary electrons we derive their spectrum from the 
gamma-ray data from Sgr B2 whose mass is 2 x 10® Mq. The derived angular distribution 


of the 330 MHz emission for 
taken from 


Jones et al. 


( 120111 ). 


his case is shown in Figure [5] together with experimental data 


We see that even for the most favourable parameters of the model the nonthermal 
component is below the intensity of thermal emission. 

The flux of nonthermal component can be subtracted from the total radio emission 
from Sgr B2 by observation of polarization. A flux of synchrotron emission is pola rized if 


the strength of large scale magnetic helds is high enough (see, e.g., 


Ginzburg 


1989|). 


However as we concluded in Section [3] the chaotic component of the magnetic held, 6H, 
is hve times larger than the large scale magnetic held. Ho, whose st rength is estimated as 


OOuG. The expected linear polarization degree, pi, was derived by 


Korchak fc Syrovatskii 


(1961), which is 

(7 + 1) (7 + 3)(7 + 5) 


Pi = 


(7 + r/3) 


32 


1 - 


( 7 ^ + 87 + 3 ) Hi 


24 


6H^ 


HI 

5H^ 


for 5H^H 


0 • 


(47) 


For the expected spectral index of secondary component 7 (for the electron spectrum 
oc E~'^) with 7 ~ 3 — 4, the estimated polarization degree of the nonthermal component of 
Sgr B2 is about pi ~ 3.5 — 4.8%. 

If, however, the turbulence inside the clouds is not strong enough and the amplitude of 
magnetic huctuations is small {SH Hq), then the polarization degree is described by the 
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equation fsee lKorchak fc Svrovatskiilll96ll) 


Pi = 


(7 + 1) 


1 - 


2 6H‘^ 


3 Hi 


for 6H < H 


0 • 


(7 + 7/3) 

The upper limit of linear polarization of synchrotron radiation can reach the level about 
75 — 80%, if the scale of Hq is about the cloud radius. 


(48) 


We hope t 


rat t 


Dickinson et al 


2ni5h 


le next generation radio telescope SKA (see IStrong et ahl 12014 
may be able to detect the nonthermal flux of Sgr B2. 


7. Hydrogen ionization by primary electrons in Sgr B2 


see 


Alternatively ionization of the molecular gas may be caused by high energy electrons 


Yusef-Zadeh et ah 


20131) . Below we analyse whether the molecular hydrogen in Sgr B2 


is ionized by electrons and whether we can distinguish from observations between processes 
of proton and electron ionization there. 


As follows from Figure [Usubrelativistic CRs are absorbed in the outer envelope, and 
electrons with energies above 1 MeV are able to penetrate into the dense core. Therefore, 
for estimates of processes of hydrogen ionization and bremsstrahlung emission from there 
we can use equations in the relativistic energy range, which simplihes the calculations 
signihcantly. The intensity of bremsstrahlung emission is determined by a cutoff position 
derived from the rate of ionization losses in the outer envelope (see, e.g.. Figure [2]). 


Electrons lose about 42 eV per one act of ionization (see 
can estimate the ionization rate of hydrogen as 

^ Up dE 

c = 


Dalgarno et al 


1999 1. We we 


42 eV dt ’ 


(49) 


where Ug is the density of relativistic electrons in Sgr B2, dE/dt represents ionization losses 
of these electrons per H-atom and per one electron. We notice that in the relativistic 
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energy range the rate of ionization losses weakly depends on the electron energy. From the 
average ionization rate in Sgr B2, ('q = 4 x 10“^® s“^, we obtain that the required density of 
electrons with any energy above 1 MeV is — 7.7 x 10“® cm“^. 

These electrons provide X-ray bremsstrahlung radiation. However, unlike the case 
of subrelativistic protons, the energy of bremsstrahlung photons is . The X-ray 


Hand ( 19981 1 


emissivity of relativistic electrons with Eg > 1 MeV can be estimated from 
For the range of 20-60 keV the bremsstrahlung emissivity Cep = 5 x 10“^^ erg 
per electron for scattering by background protons and Cee = 8.5 x 10“^^ erg for 

scattering on background electrons. The total 20-60 keV X-ray flux from Sgr B2 is 


^ -|- 6ee)Vf/f ^ 


1.7 X 10 erg s ^ cm ^ x 


(. 


M, 


H 


Vio^Mq 


(50) 


where Mh is the total mass of Sgr B2. 


For the mass of Sgr B2 M = 8 x 10^ Mfr, this X-ray flux reaches the level of 1 mCrab 


Just this flux was observed by flTerrier et al 


2010 l ). As we noticed above the origin of 


time variable continuum X-ray emission is the irradiation of Sgr B2 by hard X-ray emitted 


i n the past by S 


gr A*. As follows 


(INobnkawa et al. 


2011 


Terrier et al 


rom 


he observations of Suzaku and INTEGRAL 
2010l ) the X-ray flux is a power-law with E^ oc 


The electron bremsstrahlung of relativistic electrons is hard, oc E~^, as it follows from 
Equations fj43ll & fj4^ . Then the effect of electron bremsstrahlung is evident in the total 
X-ray flux from Sgr B2 at relatively high values of E^ as a stationary excess above see 
Figure [6l We can conclude from our calculations that the origin of ionization in Sgr B2 can 
be, in principle, determined from the INTEGRAL observations. If these observations show 
an excess of hard X-ray emission above the power-law oc at large energies, then it is in 
favour of the electron origin of hydrogen ionization. If not, then only the proton origin is 
acceptable. 
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8. Conclusion 

The conclusion of this paper can be itemized as follows: 

• The total mass and the gas distribution in the cloud Sgr B2 is highly uncertain. 
Therefore we study two extreme cases: (1) when the mass is mainly distributed in 
the relative low density outer envelope ~ 10^ — lO'^ cm“^) with the radius about 
22 pc (Model 1), (2) when the most of the mass is concentrated in the dense core 
{riH 10^ cm with the radius about 5 pc (Model 11). 

• We conclude that propagation of high energy particle in dense molecular clouds is 
determined by the amplitude of magnetic fluctuations generated by turbulent motions 
of neutral gas and their damping due to ion-neutral friction. In the outer envelope of 
Sgr B2 (if it is) with the density about 10^ — 10^ cm“^ the kinetic energy of turbulent 
motions is quite low and it does not generate fluctuations with amplitude higher than 
the strength of large scale magnetic held, Hq ~ 100 /iG. Therefore we assume that 
particles propagate their without scattering. 

• In the dense core with the density about 10® cm“® the neutral gas turbulence generates 
huctuations with the amplitude 6H ^ Hq. Our estimations show that 6H ~ 500 yuG 
for the parameters of Sgr B2 core. Then, unlike previous models of GR propagation 
in molecular clouds, we conclude that this process is determined by chaotic magnetic 
helds generated by turbulent motions of the molecular gas in the core. The key 
parameter of this process is the correlation length of the chaotic held, Lcorr which 
determines dihusion propagation of particles inside the core. For the magnetized 
particles the dihusion coefficient is constant in the relativistic energy range, and is 
proportional to y/E in the subrelativistic energy range. For the parameters of Sgr B2 
the correlation length is about Tcorr ~ 0.4 pc. 
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• For the model of hydrogen ionization by subrelativistic protons we derived their 
spectrnm depending on the hydrogen density in the onter envelope. We showed that 
the hydrogen ionization is provided by protons with the energy above 20 MeV. Their 
energy density is abont or above 50 eV cm“^. 

• The spectrnm of relativistic protons was derived from the observed gamma-ray flux of 
Sgr B2. If the mass of the cloud is mainly distributed in the outer envelope with the 
radius of 22 pc, then the energy density of these protons is about the same as near 
Earth, i.e., 1—4 eV cm“^. If, however, most of mass is concentrated in the core with 
the radius 5 pc then the required energy density of the protons is about 10 eV cm“^. 

• We conclude that if subrelativistic protons are responsible for hydrogen ionization in 
Sgr B2, they generate also the 6.4 keV line emission and continuum emission in the 
X-ray energy range < 30 keV. The expected flux in the K-a iron line from Sgr B2 is 
about 10“^ ph cm“^s“^, which is only three times below the line flux observed in 2013 
by Suzaku. It was shown that the observed flux of the line from Sgr B2 was generated 
by X-ray photons emitted by Sgr A* in the past. However, this flux is decreasing 
rapidly and we expect that in near future we might observe a component of the 6.4 
keV emission from Sgr B2 generated by CRs. 

• The expected continuum X-ray emission generated by proton bremsstrahlung is about 

= (4 — 6) X 10“^^ erg cm“^ s“^ in the range 2 — 10 keV, and = 7.5 x 10“^^ 

erg cm“^ s“^ in the range 20 — 60 keV. Both fluxes are below the continuum fluxes 
from Sgr B2 measured by Suzaku and INTEGRAL in 2009. However we expect that 
the continuum X-ray fluxes might decrease to the level predicted by the model of 
proton bremsstrahlung. It is highly desirable to get updated results from Suzaku and 
INTEGRAL. 


Subrelativistic protons generate also de-excitation gamma-ray lines. However, the 
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estimated flux of the C and O lines from Sgr B2 is too low for detection by present 
gamma-ray telescopes. We hope that it may be a target for future gamma-ray 
missions. 

• Relativistic protons in Sgr B2 produce secondary electrons by pp collisions whose 
density can be derived from the observed gamma ray flux. Thus, a flux on nonthermal 
synchrotron radio emission is expected from Sgr B2. However, observations showed 
that the flux of thermal radio emission was much higher than the estimated flux 

of synchrotron emission. Therefore, it is problematic to subtract the nonthermal 
component from the total Sgr B2 radio flux. However, unlike the thermal component, 
the flux of synchrotron emission is polarized. We estimated the expected degree of 
polarization of the nonthermal component. If the amplitude of the chaotic magnetic 
fields excited by the turbulence of neutral gas is higher than the strength of large 
scale magnetic fields in the halo, the degree of polarization is small, around 3 — 4%. 
If the kinetic energy of turbulent motions is relatively small, the large scale fields are 
not perturbed significantly by the turbulence. In this case the upper limit of linear 
polarization of the synchrotron component may reach the level of 75%. 

• An alternative model for hydrogen ionization in Sgr B2 is irradiation by high energy 
electrons. The energy of these electrons should be around or above 1 MeV. These 
electrons generate a bremsstrahlung X-ray emission with a spectrum oc E~^ which is 
harder than the emission from the Compton echo oc E~^ as observed by Suzakn and 
INTEGRAL. Therefore, we expect that the bremsstrahlung component can be seen at 
high energies. Our estimations show that the calculated bremsstrahlung component 
should be seen even in the 2009 data in the range 20-60 keV if it is generated by 
electrons. If observations do not show any excess above the spectrum E~‘^, then it 
is an argument against the leptonic origin of hydrogen ionization at least in Sgr B2, 


although more detailed analysis of the X-ray data after 2009 is needed. 
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Fig. 4.— Spectrum of subrelativistic protons in the intercloud medium, needed for the rate 
of ionization inside and outside Sgr B2 for different values of F'max- The number in the 
brackets represent the estimated energy density of subrelativistic protons in eV cm“^. The 
solid line (left side of the hgure) represents spectra of subrelativistic protons for hard spectra 
with 6 = 0.2 that gives the value of F'max = 0.2 GeV (see Fig. [2]). In this case almost all 
particles freely penetrate through the envelope, and therefore the the required spectrum of 
subrelativistic protons is the same independently of uh- The other three curve are shown 
for F'max = 0.5 GeV and different values of uh- As follow from Figure |2]for this value of 
-Fmax the spectra are relatively soft and therefore a part of particle is absorbed in the outer 
envelope. This makes the result depended on hh as shown in the hgure. The spectra of 
relativistic protons are shown by shaded area for the case of outer envelope and the line 
when the mass is concentrated in the core. Gomments are presented in the text. 
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Fig. 5.— Upper limit on the angular distr ibution of t 
line) together with experimental data from I, Tones et al. 
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ig. 6.— Spectrum of hard X-ray emission observed by INTEGRAL from 


Terrier et al. 


2010h . The bremsstrahlung emission of 10 MeV electrons for upper 1C 


Mq limits of the mass of Sgr B2 (solid lines). The X-ray spectrum form 
is shown by dashed lines. 
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